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Abstract 

The magnetized Iron CALorimeter detector (ICAL) which is proposed to be built 
in the India-based Neutrino Observatory (INO) laboratory, aims to study atmospheric 
neutrino oscillations primarily through charged current interactions of muon neutrinos 
and anti-neutrinos with the detector. The response of muons and charge identification 
efficiency, angle and energy resolution as a function of muon momentum and direction 
are studied from GEANT4-based simulations in the peripheral regions of the detector. 
This completes the characterisation of ICAL with respect to muons over the entire 
detector and has implications for the sensitivity of ICAL to the oscillation parameters 
and mass hierarchy compared to the studies where only the resolutions and efficiencies 
of the central region of ICAL were assumed for the entire detector. 

Selection criteria for track reconstruction in the peripheral region of the detector were 
determined from the detector response. On applying these, for the 1-20 CeV energy 
region of interest for mass hierarchy studies, an average angle-dependent momentum 
resolution of 15-24%, reconstruction efficiency of about 60-70% and a correct charge 
identification of about 97% of the reconstructed muons were obtained. In addition, muon 
response at higher energies upto 50 CeV was studied as relevant for understanding the 
response to so-called rock muons and cosmic ray muons. An angular resolution of better 
than a degree for muon energies greater than 4 CeV was obtained in the peripheral 
regions, which is the same as that in the central region. 


*email; kanishka@puhep.res.iii 
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1 Introduction 


The India-based Neutrino Observatory (INO) |T] is the proposed underground facility that will 
house a magnetized Iron CALorimeter detector (ICAL) designed to study neutrino oscillations 
with atmospheric muon neutrinos. In particular, ICAL will focus on measuring precisely the 
neutrino oscillation parameters including the sign of the 2-3 mass-squared difference Am ‘^2 
(= — m 2 ) and hence help to determine the neutrino mass hierarchy through matter effects. 

Oscillation signatures for neutrinos and anti-neutrinos are different in the presence of matter 
effects which become important in the few GeV energy region. These parameters are sensitive 
to the momentum P and the zenith angle cos 6 (through path length travelled) of neutrinos. 
Reconstruction of these parameters further depends on the energy and direction of muons 
and hadrons [2] produced in charge-current interactions of the neutrinos in the detector; 
hence studies of muon energy and direction resolutions are crucial. 

Since ICAL is designed to be mostly sensitive to muons, the main physics issues that ICAL 
will probe will be through studies of charged current (CC) muon neutrino (anti-neutrino) 
interactions with muons (anti-muons) produced in the final state. Two types of interactions 
are relevant: in the first, the neutrino enters the detector, interacting with (dominantly) 
nucleons in the iron nucleus. These events are identified by vertices which are inside the 
detector (tracks begin inside the detector). In the second type of interaction, the neutrino 
interacts with rock around the detector and the produced muons lose energy while propagating 
through the rock (these are the so-called rock muons or upward-going muons). These events 
have vertices outside the detector with their tracks starting at an edge of the detector. The 
first type of interaction dominantly gives low energy (few GeV) muons due to the rapidly 
falling atmospheric neutrino flux. While in the case of rock muons, most of the lower energy 
muons are stopped in the rock itself, so that the fraction of higher energy (> 10 GeV) muons 
is higher in this sample. In addition, cosmic ray muons also enter the detector from above. 

In an earlier paper [3], the response of IGAL to few-GeV muons with respect to both 
momentum magnitude and direction reconstruction was studied through simulations for muons 
generated in the central region of IGAL where the magnetic field is largely uniform both in 
direction and magnitude. Here we present for the first time a GEANT4-based simulations 
study of the muon response in the peripheral region of IGAL, where the magnetic field is 
not only highly non-uniform in both magnitude and direction, but there are signihcant edge 
effects as most of the tracks will be partially contained. Note that a substantial fraction of 
rock muon events traverses such regions in the detector; hence it is important to understand 
the muon response in these regions for such physics studies. 

The inclusion of muon response in the peripheral region of the detector can significantly 
change the oscillation and mass hierarchy results since it contains resolutions and efficiencies 
in the energy range 1-50 GeV. The first set of simulations results for precision measurement of 
neutrino oscillation parameters and the mass hierarchy were produced using only the central 
region resolutions of about 9-14% and efficiency of about 80% (see Ref. [1], [5] for more details) 
in the energy range 1-20 GeV. 

The paper is organised as follows: in Section 2 we briefly discuss some relevant details 
about the GEANT4-based simulation of the detector geometry and magnetic held, as well 
as the methodology of hit and track generation and track reconstruction. In Section 3, we 
discuss the general features of muon propagation in the different regions of IGAL. In Section 
4, we discuss the selection criteria used. We then present in Section 5 the results, with these 
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selection criteria, for the muon efficiencies and resolutions in the peripheral region of ICAL. 
A comparison of the response in all the regions of ICAL (central and peripheral) is given in 
Section 6. We conclude in Section 7, with discussions. 


2 The ICAL Detector Simulation Framework 

Details of the ICAL detector can be found elsewhere [5]. Here we briefly review the relevant 
simulation inputs for the geometry and magnetic held. The ICAL detector geometry is simu¬ 
lated [3] using GEANT4 software [6]. It comprises of three identical modules of dimension 16 
m X 16 m X 14.45 m. The direction along which the modules are placed is chosen to be the 
x-direction (and is labelled with the azimuthal angle 0 = 0) and the direction perpendicular 
to the x-axis in the horizontal plane is the ^/-direction. The vertical direction is the ^-direction 
with the z-axis pointing upwards (so the polar angle is also the zenith angle). Each module 
comprises a stack of 151 layers of 5.6 cm thick magnetized iron, separated by a 4 cm air gap 
in which the active detector elements, the RPCs are placed. The |/-direction is in the plane 
of the iron plates, parallel to the slots of the magnet coil, as shown in Fig. [H with the origin 
of the coordinate system being the centre of the central module. 

Apart from coil slots (thin 8 m long slots centred around ?/ = 0 at x = a;o±4 m, where Xq is 
the central x-coordinate of each module), there are vertical steel support structures which are 
placed at every 2 m along the whole detector in both x and y directions, to maintain the air 
gap between plates. These are dead spaces that affect the muon reconstruction. In addition, 
the magnetic held is also not uniform everywhere and so the quality of reconstruction depends 
on the region where the event is located. 

2.1 The Magnetic Field 

The magnetic held has been simulated in a single iron plate using the MAGNET6 software 
[7]. The magnetic held map is generated at the centre (z = 0) of the plate and is assumed 
uniform over the entire thickness of the plate. The magnetic held lines in a single iron plate 
in the central module are shown in Fig. [1] The held is generated by current circulating in 
copper coils that pass through slots in the magnetized plates. The slots can be seen as thin 
white vertical lines at x = ±4 m in the hgure. The direction and length of arrows denote the 
direction and magnitude of the magnetic held. 

The “central region” is dehned as the volume contained within the region |x| < 4 m, \y\ <4 
m with 2 ; unconstrained, that is, the region within the coils slots in each module. The central 
region has the highest magnetic held that is uniform in both magnitude and direction {By) 
while the region labelled “peripheral region” (outside the central region in the y direction, 
\y\ > 4 m) has the maximally changing magnetic held in both direction and magnitude, with 
the held falling to zero at the corners of the module. 

In an earlier paper [3] , the muon response was studied in the central region; here we study 
the peripheral region where, apart from the changing magnetic held, the reconstruction is also 
ahected by edge ehects. In addition, there are two small regions outside the coil slots in the x 
direction (|x| >4 m), labelled as the “side region” in Fig. [1] where the magnetic held is about 
15% smaller than in the central region and in the opposite direction. The side regions in the 
central module are contiguous with the side regions in the adjacent modules and the quality 
of reconstruction here is expected to be similar to that in the central region. However the left 
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Fig. 1: Magnetic field map as generated by the MAGNET6 [7j software in the central iron layer 
of the central modnle. Points A = (0, —650, 0) cm, B = (300, —650, 0) cm, C = (—2270, 0, 0) 
cm (in the 1^* modnle of ICAL), are dehned for later use. Notice that C is actually in the 
left-most module of the detector and is simply marked here for convenience. 


(right) side region of the left (right)-most module will suffer from edge effects and we shall 
consider them separately in the study. 

2.2 Event Reconstruction 

In each analysis, 10,000 muons are propagated in the detector with hxed momenta and di¬ 
rection and with the starting point of the muons in either the peripheral or side regions. In 
contrast to the earlier study in the central region [3], here the muon momenta vary from 1-50 
GeV/c, with the higher energies being of interest for rock muon studies. The muon tracks are 
reconstructed based on a Kalman hlter [8] algorithm. The RPGs that signal the passage of 
muons through them have a position sensitivity of about ±1 cm in the x- and ^/-directions (the 
RPG strip width is 1.96 cm) and about ±1 mm in the ^-direction (the gas gap in the RPGs 
being 2 mm). In addition, “hits” or signals can be generated in adjoining RPG strips as well, 
and further that the RPG efficiency and time resolution are about 95 % and 1 ns respectively 
[9]. For more details regarding the generation of hits, tracks and their reconstruction, see 
Ref. [3]. 

3 General Feature of Muon Response in the Peripheral 
and Side Regions 

We hrst discuss the general expectations, based on the detector geometry and the orientation 
of the magnetic held. 

The Lorentz force equations are F = q{v x B), where B is the magnetic held that is 
conhned to the x-y plane, q is the charge of the particle with momentum P and energy E 
so that its velocity v is directed along the momentum with magnitude v = Pc^/E. Since 
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g = — 1 for /i , the components of force in the peripheral region for an upward-going muon, 
momentarily ignoring energy loss, are. 


Fx 



Fy 

"^zBx 1 


F. 

^yBx '^xBy , 

(1) 


whereas the analogous components of force in the side region are given as: 

Fx V^By , 

= 0 ; 

F^ = VxBy . (2) 

It is seen that in both regions, Fx and Fy are independent of 0 (that is, independent of the 
momentum components in the plane of the held) and depend on (i.e., on the zenith angle 
9 alone) while F^ depends on both 9 and 0. Depending on the components of magnetic held 
Bx and By, Eqs. [T] and [2] give the net force in the diherent regions of ICAL. Consider the 
in-plane (x, y components) forces in the regions denoted as 1-10 in Fig. [2l It can be seen that 
in regions 1, 2, 7, 8, Fy is such that it causes an upward-going muon (cos6' > 0) to be bent in 
a direction going out of the detector, thus lowering the reconstruction efficiency. The ehect 
is just the opposite in regions 3,4, 5, 6. If F^ > 0 as well, the already upward-going muon 
traverses more iron layers as discussed in Ref. [3] and hence gives good resolution; hence, Fy 
ahects the reconstruction efficiency while F^ determines the quality of reconstruction. Since 
Fz depends on both {9, 0) as well as the magnetic held, the sign of F^ is shown in Fig. |2] 
inside a circle of 0 in each of regions 1,2, 3,4 (for negatively charged upward-going muons), 
for I Fa; I ~ \By\ with purple (cyan) regions denoting F^ > (<)0. 

Hence the magnetic held that determines the quality of reconstruction breaks the azimuthal 
symmetry, so muons in diherent 0 directions (for the same momenta and polar angle cos 6*) 
have diherent detector response. This was discussed in detail in Ref. [3]. Going by the 
force equations, we therefore analyse the muon response in the peripheral region in four 
diherent set of 0 bins as shown in Fig. |31 bin I: |0| < 7r/4, bin II; 7r/4 < 0 < 37r/4, bin III; 
—37r/4 < 0 < —7r/4, and bin IV: 37r/4 < |0| < vr. 

For muons with starting point in the negative y peripheral region (regions marked 1, 2, 3,4 
in Fig. [2]), this implies that most of the muons with momenta such that 0 lies in bin III (but 
otherwise having same magnitude and cos6') are prone to exit the detector from the side; 
however, there will be marked diherences in the quality of reconstruction between regions 
(1,2) and (3,4) due to the diherent Fy force that turns the track back into the detector in 
regions 3,4. Hence the average detector response in this region is an average over these two 
diherent behaviours. In addition, F^ > 0 (< 0) for bin H muons in regions 1,2 (3,4) and this 
helps to improve the reconstruction in regions 1,2, so that bin H muons can be expected to 
have the best quality of reconstruction of the regions 1-4. 

A similar analysis can be done for muons which begin in the positive y peripheral region 
(regions marked 5, 6, 7, 8). Of course, tracks at the edge of a region or of high energy muons 
may move from one region to another where a diherent magnetic held applies; however, we 
simply bin the events according to the region in which the muon originates. 
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Fig. 2: Magnetic field map with the net force directions in the peripheral and side region. The 
thick black arrows indicate the direction of the magnetic held with labels , 

that denote the sign of the B^, By components in each region. The brown arrows indicate the 
direction of or Fy force components that will act on a negatively charged upward-going 
muon. The small coloured circles indicate the direction of F^ in each region, with purple 
(cyan) denoting F^ > 0(< 0). (Note that side regions 9 and 10 are in the and 3^'^ modules 
of the detector respectively and are shown together in the same module for convenience.) 



X 



Fig. 3: The choice of 0 bins in the peripheral (left) and side (right) regions. 
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In the side region 9, causes the particle to go out of the detector but > 0 and so 
this region has good resolution but worse efficiency. The results are opposite in side region 
10 since By < 0. We therefore dehne the same (j) bins for the side region as were used for 
the central region, viz., bin I: \(j)\ < vr/d, bin II: vr/d < |0| < 7r/2, bin III: 7r/2 < \(j)\ < Svr/d, 
and bin IV: 37r/d < |0| < vr. The difference is in the dehnition of the second and third bins, 
see Fig. [3l and is more appropriate from the point of view of the geometrical conhguration 
in this region. Region 9 (10) will have the worst reconstruction in 0 bin IV (I). However, in 
region 10, the direction of the force is expected to improve the results just as in the case 
of regions 3, d. 

The results will be the same for downward-coming (with cos 6* < 0) while that for 

downward-coming muons or upward-going anti-muons can be obtained by symmetry. In our 
analysis, therefore, we study the muon response in the peripheral regions 1-d and the side 
region 9. 

It is important to keep in mind that the support structures and coil gaps also break this 
azimuthal symmetry in a non-trivial way and the effects of the geometry may alter the trends 
of the distributions as discussed above. 

The net effect of the nature of the detector geometry and magnetic held can be seen in 
the muon resolutions and efficiencies and we will now discuss these. 

As discussed above, we study the response in the following peripheral and side regions. 

In the Peripheral Region : Here, 10,000 muons (/i“) were propagated with hxed input 
momenta Pin and direction cos6' (and smeared over the entire azimuthal angle 0), with their 
starting point uniformly smeared over the region centred at (0, —600, 0) cm and extending 
upto ± (2400, 200, 720) cm from it; this comprises the whole peripheral region along the three 
modules of the detector in the negative y region where the magnetic held is non-uniform. 

In the Side Region : In the side region, muons (/i“) were propagated with the same 
procedure as above but with point of origin smeared in a region centred around (—2200, 0, 0) 
cm and (2200, 0, 0) cm (which are in the R* and 3^'^ modules of the detector respectively) and 
smeared uniformly in ± (200,400, 720) cm around these. 


4 Selection Criteria Used 

All tracks which satisfy the loose selection criterion y^/ndf < 10, are used in the analysis, 
where is fhe standard y-squared of the ht and ndf are the number of degrees of freedom, 
ndf = 2 X Nhits — 5, where Nhus are the number of hits in the event, Nhus > 5 and the Kalman 
hlter involves the htting of 5 parameters. 

Further selection criteria are used to get reasonable hts and hence resolutions. Two major 
constraints have been applied in both the peripheral and side regions to remove low energy 
tails. The hrst is similar to that applied when analysing tracks in the central region [3]: the 
Kalman hlter algorithm may generate more than one track. While this may be correct and 
useful in the case of genuine neutrino CC interactions, where one or more hadrons accompany 
the muon, this is a problem for single muon analysis and arises because of detector dead 
spaces (for instance, two portions of a track on either side of a support structure may be 
reconstructed as two different tracks). This problem will be mitigated by the identihcation 
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of a vertex in a genuine neutrino interaction, here, we place a constraint and analyse only 
those events for which exactly one track is reconstructed, leading to a consequent loss in 
reconstruction efficiency. The second selection criterion is specihc to the peripheral and side 
regions and is described below. 

Initially, events were generated at hxed points of origin to understand the effect of the mag¬ 
netic held. In the peripheral region, the starting point was chosen to be either at point A (in a 
region of nearly zero magnetic held) or B (large magnetic held with both x- and ^/-components 
non-zero), while in the side region, a generic point C was studied (see Fig. [T]). Results of this 
study clearly indicated that a large fraction of events whose tracks were truncated because 
the particle exited the detector (so-called partially contained events) were relatively poorly 
reconstructed. These could not be eliminated by tightening the constraint on of fho hts; 
however, they could largely be removed by demanding that the track contain a minimum 
number of hits, such that either > uq or Nhus/ cos 9 > uq (note that there may be 

multiple hits per layer), where Uq needed to be carefully optimised. It was found that for a 
given momentum and direction of the muon, hq needed to be larger (smaller) in regions where 
the magnetic held strength is small (large). Where the muon does not leave the detector, so 
that the entire track is contained in the detector (fully contained events), no constraint on 
Nhitg is needed. With this understanding, the generic peripheral and side region response was 
studied. 

4.1 Effect of Selection Criteria 

The effect of Nhug > uq or Nhus/cos 9 > no can be seen from Fig. IHfor the peripheral region. 
If the event is fully contained, there is no constraint on Nhug] the effect of Uq = 15 is shown 
in the left-hand side of Fig. 0] for (Pj„,cos6*) = (5 GeV/c, 0.65) and (9 GeV/c,0.85) where 
the histogram in the magnitude of the reconstructed momentum P^ec is plotted. For Pin = 5 
GeV/c, it is noticed that the Nhus constraint does not affect the Prec momentum distribution 
much, as most of the events are fully contained. But it gives a better (more symmetrical) shape 
to the distribution by reducing the low-energy tail. On the other hand, the effect for Pj„ = 9 
GeV/c is stronger, with the hump at lower Prec being eliminated with the Wits selection 
criteria. Fig. [5] shows the effect of Wits > no or Nhus/cos 9 > no on Wits distributions in the 
peripheral region. In all cases, the few events surviving below the constraint are from totally 
contained events, on which no constraint is placed. These cause the histograms to remain 
non-zero in the region Wits < no or Wits/ cos^ < rto as seen in Fig. 01 However, these events 
are relatively few in number, being less than 2% (3%) of the total reconstructed events for 
no = 15 (20) in Fig. [5l The constraint Nhus/ cos 9 > 15 is the most conservative one, with 
a loss of only about 10% of the reconstructed events and is found to be more optimal than 
Nhits ^ no- 

Different choices of no can be used. We have shown the effect of (a) no constraint, (b) 
Wits > 15, and (c) Wits/ cos9 > 15 in the left-hand side of Fig. 01 The last choice is motivated 
by the fact that a slant-moving muon (in the absence of magnetic held) would move a distance 
d/ cos 9 in comparison to a vertically upward-going muon of the same momentum that would 
traverse a distance d. Similar hgures on the right use the choice no = 20, with the more 
stringent requirement giving distributions with correspondingly smaller root-mean-square or 
square root of the variance (RMS widths) by about 7-8%, but showing a decrease in the total 
number of reconstructed events by 10-15%. Note also that increasing no eventually leads to 






Fig. 4: Top (bottom) figures show the reconstructed momenta P^ec using selection criteria 
Nhits > uo for partially contained events in the peripheral region with {Pin, cos 9) = (5 
GeV/c, 0.65) (top) and (9 GeV/c, 0.85) (bottom) with uq = 15 (20) in the left (right) figure. 
Fully contained events have no Nhns constraint. In each figure, the black curve is without 
constraints on Nhus, red is with Nhus/cos 9 > Uq and blue is for N^us > riQ. 
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Fig. 5: Top (bottom) figures show the Nhus distributions using selection criteria Nhus > for 
partially contained events in the peripheral region with (Pj^, cos6') = (5 GeV/c, 0.65) (top) 
and (9 GeV/c, 0.85) (bottom) with no = 15 (20) in the left (right) figure. Fully contained 
events have no Nhns constraint. In each figure, the black curve is without constraints on Nh^s, 
red is with Nhus/ cos 9 > uq and blue is for N^its > ^q. 
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removal of well-reconstructed events, as visible from the loss of events in the peak apart from 
just trimming the tails for the lower momentum = 5 GeV/c for hq = 20. 

Similarly, the effect of the selection criteria on the reconstruction in the side regions is 
shown in Fig. O Here the constraint on the partially contained events is not as strongly 
marked as in the peripheral region: while there is certainly a decrease in the RMS width of 
the distribution and in the number of selected events when the constraint is applied, a larger 
fraction of events are lost due to the constraint, with a number of “good” events being lost 
from the peak of the distribution as well, unlike in the peripheral region. 




Fig. 6: The hgures show the reconstructed momenta Pj-ec using the selection criteria 
for partially contained events in the side regions 9 (left) and 10 (right) for (Pin, cos0) = (5 
GeV/c, 0.65) with no = 15. Fully contained events have no Nhus constraint. In each hgure, 
the black curve is without constraints on Nhus, red is with Nhus/ cos 9 > uq and blue is for 
^hits ^ rig- 


The hnal choice of selection criteria will be guided by the physics study. In case the 
requirement is good momentum resolution, then the choice uq = 20 may be appropriate (that 
is, either Nhus > 20 or Nhus/ cos 9 > 20). However, since the shape of the distribution is 
already reasonable for rig = 15, this choice may be used when the focus is not so much on 
precision reconstruction but on higher event reconstruction rates. In the rest of this paper, we 
shall apply the constraint Nhus/ cos 9 > 15 as being appropriate and sufficient. This choice 
also improves the reconstruction efficiency of large angle (small cos 6*) events whose tracks 
naturally contain fewer hits and are harder to reconstruct. 

In the next section, we present the results on muon resolution and efficiencies in the 
peripheral and side region using these selection criteria. 


5 Muon Response in the Peripheral and Side Regions 

5.1 Momentum Reconstruction Efficiency 

The reconstruction efficiency is dehned as the ratio of the number of reconstructed events n^ec 
(irrespective of charge) to the total number of events, Ntotal- We have 

erec = 1^, (3) 

total 

with error = Verec(l - erec)/Wotal • 
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Fig. [7] shows the reconstruction efficiency averaged over as a function of input momentum 
for different cos 6* values in the peripheral and side regions. 




Fig. 7: Reconstruction efficiency averaged over all cj) bins as a function of the input momentum 
Pin (GeV/c) for different zenith angles cos 6* in the peripheral (left) and side 9 (right) regions. 
For a discussion of the selection criteria see the text. 

The reconstruction efficiency increases for all angles, from Pin = 1 GeV/c since the number 
of hits increases as the particle crosses more layers. Since there are fewer hits for more slant¬ 
angled muons, the efficiency at a given momentum is better for larger values of cos 6. Also, 
the reconstruction efficiency is very similar for all the peripheral and side regions. 

In all cases, the slight worsening of the efficiency for cos 6 = 0.85 at higher momenta is 
spurious and is due to the selection criterion that the event should reconstruct exactly one 
track. At large angles, it is more likely that two portions of a track on either side of a dead 
space such as a support structure are reconstructed as two separate tracks. Efforts are on to 
retrieve such events by improving the reconstruction code [lO] . When these tracks are correctly 
reconstructed, the efficiency is expected to saturate rather than fall off at these momentum 
values. Again, such tracks are not expected to be troublesome in genuine neutrino events, as 
discussed earlier. 


5.2 Relative Charge Identification Efficiency 

The charge identihcation (cid) of the particle is critical in many studies since it distinguishes 
events initiated by neutrinos and anti-neutrinos; these have different matter effects as they 
propagate through the Earth and hence give the required sensitivity to the neutrino mass 
hierarchy. The charge of the particle is determined from the direction of curvature of the track 
in the magnetic held. Relative charge identihcation efficiency is dehned as the ratio of the 
number of events with correct charge identihcation, ndd, to the total number of reconstructed 
events, Urec, i.e., 

^cid 

*^cid — 

^rec 

where the errors in ndd and Urec are correlated so that the error in the ratio is calculated as: 

<^Ccid '\/^cid(l *^cid)/Wec • 
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Fig. [8] shows the relative charge identihcation efficiency as a function of input momentum 
for different cos 6* values in the peripheral and side region 9. (Similar results apply for side 
region 10). The muon undergoes multiple scattering while propagating in the detector; for 
small momentum, since the number of layers traversed is small, this may lead to an incorrectly 
reconstructed direction of bending, resulting in the wrong charge identihcation. Hence the 
charge identihcation efficiency is relatively poor at lower energies but as the energy increases 
cid efficiency also improves. At very high input momenta, bending due to the magnetic 
held is less. For partially contained events, only the initial relatively straight portion of the 
track is contained within the detector; this leads to large momentum uncertainty as well as 
mis-identihcation of charge. Overall the relative charge identihcation efficiency is marginally 
smaller than in the central region because of the smaller magnetic held. 




Fig. 8: Charge identihcation (cid) efficiency averaged over all 0 bins as a function of the input 
momentum (GeV/c) for diherent zenith angles cos 6 in the peripheral (left) and side 9 
(right) regions. For a discussion of the selection criteria see the text. 


5.3 Direction (up/down) Reconstruction 

The reconstructed zenith angle distributions for Pin = 1 GeV/c at cos^ = 0.35 and cos 9 = 
0.85 in the peripheral and side region 9 are shown in Figs. P and fTOl respectively. 

From Fig. [HI it is noticed that there are few events reconstructed in the downward direction 
(wrong direction) with 9rec > 7r/2. For Pin = 1 GeV/c with cos6' = 0.35 (0.85), this fraction 
is about 0.48 (0.89)% and it drops to a negligible value at higher energies for all cos6*. Similar 
results are obtained for the side region 9 as can be seen from Fig. [101 This small fraction also 
contributes to wrong cid since the relative bending in the magnetic held is measured w.r.t the 
muon momentum direction. 

The direction determination depends on the time resolution while the charge identihcation 
depends also on the strength of the magnetic held. A 1 GeV/c muon with cos 0 ~ 1 traverses 
about 12 layers; this corresponds to a time diherence between hrst and last hit of about 4 ns. 
Since the RPGs have a time resolution of 1 ns, this explains why the fraction of muons whose 
direction is wrongly determined is small. 
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Fig. 9: Reconstructed zenith angle distributions for = 1 GeV/c at cos 9 = 0.35 (left) and 
0.85 (right) respectively, in the peripheral region. Note that the y-axis scales are different for 
the two plots. 



Fig. 10: Reconstructed zenith angle distributions for Pin = 1 GeV/c at cos 6 = 0.35 (left) and 
0.85 (right) respectively, in the side region 9. Note that the y-axis scales are different for the 
two plots. 
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5.4 Zenith Angle Resolution 

Those events which are successfully reconstructed (for all 0) are analysed for their zenith 
angle resolution. The events distribution as a function of the reconstructed zenith angle 6rec 
is shown in Fig. [TT] for a sample input (Pin, cos 6^) = (5 GeV/c, 0.65) for the peripheral and 
side region 9 respectively. 




Fig. 11; Reconstructed distribution 9rec for input (Pin, cos 0) = (5 GeV/c,0.65) in the periph¬ 
eral region (left) and side region 9 (right). The selection criteria are the same as before. 

The angular resolution is good in both the regions and is in fact better than about a degree 
for input momentum greater than a few GeV, as seen in Fig. [T21 with the resolution being 
marginally better in the side region. Similar results are obtained in side region 10 as well. 
In addition, the fraction of events reconstructed in the wrong direction (wrong quadrant of 
cos 9) is negligibly small, being less than 0.5% for Pin > 2 GeV/c. 




Fig. 12: Resolution, ag, of reconstructed angle 9rec as a function of the input momentum Pjn 
(GeV/c) for different values of input cos 9 in the peripheral region (left) and side region 9 
(right). The selection criteria are the same as before. 
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5.5 Muon Momentum Response 

While the cid efficiency and zenith angle resolntion are insensitive to the azimnthal angle 
(j), dne to the reasons given above, we analyse the mnon momentnm response in different 
0 bins. The response is shown in Fig. [13] for the peripheral region with the constraint 
Nhits/ cos 0 > 15 being applied as nsnal to the partially contained events, for sample inpnt 
valnes of {Pin, cos 9) = (5 GeV/c,0.65). 

The histograms in P^ec have been fitted with Ganssian fnnctions. The width of each 
distribntion of the fonr sets differs while the mean remains similar. As expected, 0 bin III 
(with most mnons exiting the detector from the side) has the smallest nnmber of reconstrncted 
events and the worst resolntion. Bin II has the best resolntion, while bins I and IV have a 
similar response. This is in contrast to the response in the central region [3] where the 
reconstrnction efficiencies were ronghly eqnal in all 0 bins. 






Fig. 13: Ganssian fits to reconstrncted momentnm distribntions Prec (GeV/c) for mnons with 
fixed energy {Pin, cos 6) = (5 GeV/c, 0.65) in fonr different bins of azimnthal angle in the 
peripheral region. See text for details on the bins and the selection criteria nsed. 

Similar histograms are shown in Fig. [TUfor side region 9. As discussed earlier, 0 bin IV 
has both the worst reconstruction and the worst resolution, while bin I has the best ones. 
Unlike the peripheral case where the bins I and IV had similar response, here bins II and III 
are not similar because the side region is not symmetric between these two bins: muons in 
bin III are more prone to exit the detector and hence the detector response is worse in both 
efficiency and quality of reconstruction. 

The results in region 10 are similar to region 9 with interchange of bins I and IV, and bins II 
and III as can be easily understood from Fig. [2] However, overall the quality of reconstruction 
is better in region 10 by about 15% due to the nature of the forces in this region as discussed 
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earlier; see Fig. [2l We shall show results for side region 9 everywhere, as being the more 
conservative result and simply remark on similarities/differences to be expected in region 10. 






Fig. 14: Gaussian hts to reconstructed momentum distributions P^ec (GeV/c) for muons with 
hxed energy (Pm, cos 6*) = (5 GeV/c, 0.65) in four different bins of azimuthal angle in the side 
region 9. See text for details on the bins and the selection criteria used. 

Fig. [T5] shows the momentum resolution as a function of Pm in the peripheral region for 
the four 0 bins using the selection criteria Nhus/ cos 9 > tiq with uq = 15, 20, for cos6' = 0.65. 

5.6 Momentum Resolution as a Function of (0, 0) 

Gaussian hts to the reconstructed momentum distribution in these regions give the recon¬ 
structed mean and RMS width a. The momentum resolution (P) is dehned from these hts 
as. 


R = a/Pin, (5) 

with error 5R = 6a/Pin ■ 

Fig. m] shows the variation of resolution as a function of Pm from 1 to 50 GeV/c for diherent 
values of cos 6 from 0.35 to 0.85 in the diherent 0 bins of the peripheral region. In all bins, 
the momentum resolution improves with the increase of energy upto about Pm ~ 6 GeV/c as 
the number of hits increases, but worsens at higher momenta since the particle then begins 
to exit the detector. This ehect is considerable in the 0 bin III which therefore has the worst 
resolution while 0 bin II has the best resolution, as expected from the earlier discussions. In 
general, the resolution improves for more vertical angles (larger cos 9) as the number of hits 
in a track increases. 
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Fig. 15: Muon resolution in the peripheral region as a function of input momentum Pin 
(GeV/c) for cos 9 = 0.65 in different bins of 0 with Nhus/ cos 9 > uq cut, where uq = 15 (20) 
left (right). 


Fig. [T7I shows similar results for the side region 9. Again, it is observed that for all the 
angles and energies, 0 bin I has the best response while the resolutions worsens in bins III 
and IV. Results in region 10 are similar to those in region 9 with interchange of response in 
0 bins (I, IV) and (II, III), but with a few percent better resolution in all cases. 

The resolution for a given Pin is marginally better in the side region than in the peripheral 
region due to the somewhat larger and uniform magnetic field. A detailed comparison of the 
response in different regions will be presented in the next section. 

6 Comparison of Muon Response in Different Regions 
of ICAL 

We compare the muon response in the peripheral and side regions with that in the central 
region as presented in Ref. [3]. 

For all choices of selection criteria, the reconstruction and cid efficiencies in the central 
region are better than either the peripheral or side region as shown in Fig. [T8]for cos 9 = 0.65; 
however, for input momenta upto Pj„ ~ 8 GeV/c, the central and side region cid efficiencies 
are comparable. Note that applying more stringent selection criteria in order to improve the 
momentum resolution in the peripheral and side regions (and hence overall resolution of the 
detector) will further worsen the reconstruction efficiencies in these regions. 

In addition, the angular resolutions are very similar between the peripheral and side re¬ 
gions, as can be seen from Fig. [12] and are in fact similar to those obtained earlier in the 
central region [3]. 

The comparison of the 0-averaged peripheral and side region momentum resolutions as 
a function of input momentum Pin from 1 to 15 GeV/c is shown in Fig. [T9| for cos 9 = 
0.45, 0.65, 0.85. We have also shown the 0-averaged central region results [5] in the same 
plots. The criterion of a single reconstructed track only was also applied to the central region, 
but no constraint was placed on Nhns- While the side region resolutions are only marginally 
better than those in the peripheral region, the central region gives the best resolution, as 
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Fig. 16: Muon resolution in the peripheral region as a function of input momentum Pin 
(GeV/c) for different values of cos6' in different bins of 0. 






Fig. 17: Muon resolution in the side region 9 as a function of input momentum Pin (GeV/c) 
for different values of cos 6 in different bins of <p. 
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Fig. 18: Comparison of reconstruction (left) and cid efficiency (right) of central, peripheral 
and side regions as a function of Pin (GeV/c) at cos 6* = 0.65. Note that the y-axis scales are 
different for the two plots. 


expected. However, we note that the results are 0 averaged and so the resolutions can be much 
improved in the peripheral and side regions depending on the 0 bin chosen. The peripheral 
and side region resolutions can be improved by changing the selection criteria at the cost 
of reconstruction efficiency. The resolutions in all regions are comparable at low momenta, 
-Pin < 3 GeV/c, since almost all tracks are fully contained in this case. 




2 4 6 8 10 12 14 

P,„(GeV/c) 



2 4 6 8 10 12 14 

PJGeV/c) 


Fig. 19: Gomparison of resolutions in peripheral and side region 9 as a function of the input 
momentum Pin (GeV/c) along with earlier results in the central region |3] for different values 
of COS0 = 0.45,0.65,0.85. 


7 Discussions and Conclusion 

The goal of the proposed IGAL detector is to study neutrino oscillations using atmospheric 
neutrinos. It is more sensitive to muons and hence the physics will focus on charged current 
scattering of (z7^) in the detector. Hence a simulations study of the response of IGAL 
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to muons is crucial. The ICAL geometry was simulated using GEANT4 software and the 
detector response was studied for muons with momenta from 1 to 50 GeV/c, polar angle 
cos6' > 0.35 and smeared over all azimuthal angles, —vr < (p < tt. In the current study, muons 
were generated in the peripheral and side region of the ICAL detector where the magnetic 
held is non-uniform in both magnitude and direction and where edge effects are important. 
The study showed that a crucial selection criterion on the number of hits Nhus/ cos 9 > hq for 
partially contained tracks was necessary to achieve good detection efhciency. The magnetic 
held and the detector geometry break the azimuthal symmetry; hence the muon response was 
analysed in diherent 0 bins. 

Results using Nhus/ cos 9 > 15 show that the best momentum resolutions of about 10-15% 
are obtained in bin II (vr/d < 0 < Svr/d) at input momenta of Pin > 4 GeV/c in the peripheral 
region and in bins I and II (|0| < 7r/4 and 7r/4 < |0| < 7r/2) in Side region 9 (see Figs. [2] and 
in] for dehnitions of these regions). 

Also, 0-averaged results are obtained with Nhns/ cos 9 > 15 for the reconstruction efh¬ 
ciency, charge identihcation efhciency and momentum resolution as shown in Figs. [18] and [19] 
for the peripheral and side regions of the ICAL detector in comparison with earlier results in 
the central region [3]. 

A reconstruction efhciency of about 60-70% and a correct charge identihcation of about 
97% of the reconstructed muons was obtained for Pin > 4 GeV/c and this decreased to about 
90% for higher momenta Pin ~ 50 GeV/c in both regions. Average (over 0) resolutions ob¬ 
tained are between 15-25% over Pin = 1-15 GeV/c in the peripheral region and marginally 
better in the side region, with the central region response being the best. Note that these 
responses are relevant for studies such as precision measurement of neutrino oscillation pa¬ 
rameters or the mass hierarchy determination with ICAL. For the case of physics studies such 
as rock muons or cosmic ray muons, the response in only certain 0 bins are relevant since 
the muons in these cases are always entering the detector from outside; for this reason, the 
performance will be better than the averages shown here. 

In contrast, good angular resolution of better than a degree for Pin > 4 GeV/c is obtained 
in the peripheral and side regions, which is comparable to that in the central region. 

The simulations indicate that the detector has a good response to muons, with recon¬ 
struction of momentum with 15-24% resolution, direction reconstruction of about a degree 
for muon energies greater than 4 GeV and charge identihcation of about 97%. While fully 
contained events are reconstructed with the same efhciency as in the central region, only those 
partially contained ones which have at least Nhus/cos 9 > uq in their tracks, uq ~ 15, are well 
reconstructed in the simulations. This implies a loss of reconstruction efhciency due to this 
criterion. However, the number of events reconstructed in these regions, which is expected 
to be about 50% from naive considerations of detector geometry, is about 60-70%, due to 
the ehect of the magnetic held, which increases the recontruction efhciency in the peripheral 
region. 
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